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Using a specific version of thermal Quantum Field Theory (QFT), 
supplemented by operator-field evolution of the Langevin type, 
we discuss two issues concerning the Bose Einstein correlations 
(BEC): the origin of different possible coherent behaviour of the 
emitting source and the origin of the observed shape of the BEC 
function C2(Q). We demonstrate that previous conjectures in this 
matter obtained by other approaches are confirmed and have re- 
ceived complementary explanation. 



1. Introduction 

The Bose-Einstein correlations (BEC) are since long 
time recognized as very important tool providing in- 
formation's about hadronization processes not available 
otherwise, especially in what concerns space-time exten- 
sions and coherent or chaotic character of the hadroniz- 
ing sources. Because the importance of BEC and their 
present experimental and theoretical status are widely 
known and well documented (see, for example, pp and 
references therein), we shall not repeat it here. Instead 
we shall proceed to main point of our interest here, al- 
ready mentioned above, namely to discussion of: (i) how 
the possible coherence of the hadronizing system influ- 
ences the two body BEC function C 2 (Q) 



C 2 (Q) 



N 2 {k,k') 
N x Qt)Ni{ht) 



(1) 



and (ii) what is the true origin of the experimentally 
observed Q-dependence of the C 2 (Q) correlation func- 
tion in the approach used here (out of which the space 

ha 



time information is being deduced, Q = in/Ji 



-\J (kfj, — fcp) 2 with fc M and k'^ being the four-momenta of 
detected particles; in what follows we shall for simplicity 
assume that all produced particles are bosons). In liter- 
ature one finds that in some approaches using quantum 
statistical methods [2] 



whereas in other approaches 
C 2 (Q) = l + X-Q(Q-r), 



(3) 



C 2 (Q) = 1 + 2p(l - p) ■ y/n(Q -r) +p 2 - Q(Q ■ r), (2) 



(actually, this is the most frequently used form). In 
both cases parameters p and A are called coherence pa- 
rameters defining the degree of coherence of hadroniz- 
ing source (for purely coherent source p = A = and 
there is no BEC, for p = A = 1, i.e., in purely chaot- 
ic case, both equations coincide). Although in E] 
they are operationally expressed in the same way, i.e., 
p = X = {N chao ti C ) / (N total) , we shall formally differen- 
tiate between them when using the corresponding ex- 
pressions for C 2 (Q) because, as will be shown later, 
the concepts of coherence they correspond to are dif- 
ferent in each case. The choices of Q(Q) discussed in 
literature vary between EJ © (here q = r ■ Q with 
r = = y/r^r^ being a 4- vector such that (t>) 2 has 
dimension of length and product Q ■ r = Q^r^ = q is 
dimensionless): 

• Gaussian: fl(q) = exp (— Q 2 r 2 ); 

• exponential: £!((?) = exp(— Qr); 

• Lorentzian: f2(g) = 1/ (1 + Qr) 2 ; 

• given by Bessel function p]: 

0(g) = [Ji{Qr)/{Qr)f . " ' 

The most frequently used forma are Gaussian and ex- 
ponential ones (J]. The above questions can be therefore 
rephrased in the following way: (i) why are forms of 
C 2 (Q) in eqs. (J2J and 10 different and are parameters p 
and A referring to the same quantity and (ii) what stays 
behind specific choice of the form of fl(Q) function as 
listed above. 

To make our point more clear we shall work here di- 
rectly in phase space (as, for example, in 0), no space- 
time considerations will be used (contrary to the major- 
ity of works on BEC PI El El)- ^ s our working tool we 
shall choose some specific (thermal) version of Quantum 
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Field Theory (QFT) supplemented by the operator-field 
evolution of Langevin type proposed recently [HI 03 E3 • 
We shall demonstrate that: (i) the origin of differences 
in C2(Q) in eqs. (3J and © lies in different ways of in- 
troducing the concept of coherence in both approaches, 

1. e., p and A referred to different concepts of coherence 
in each case; (u) in order to obtain a given (experimen- 
tally observed) shape of the BEC correlation function 
C2{Q) (i.e., the Q(Q)) one has to account somehow for 
the finiteness of the space-time region of the particle 
production (i.e., of the hadronizing source). In QFT ap- 
proach used here it is particularly clearly seen and is 
connected with the necessity of smearing out of some 
generalized functions (delta functions: 8{Q I1 = k ii — k'^)) 
appearing in the definition of thermal averages of some 
operators occurring here. The freedom in using different 
types of smearing functions to perform such a procedure 
allows us to account for all possible different shapes of 
hadronizing sources apparently observed by experiment. 
(Actually, careful inspection of all previous approaches 
to BEC using QFT, cf, for example, [H], shows that this 
was always the procedure used, though never expressed 
so explicitly as is done here.). 

2. Description of hadronizing source 

Let us recapitulate now the main points of our approach 
(for details see j2]). The collision process produces usu- 
ally a large number of particles out of which we select 
one (we assume for simplicity that we are dealing only 
with identical bosons) and describe it by operator b(k, t) 
(the notation is the usual one: b(k,t) is an annihilation 
operator, k is 3-momentum and t is a real time). The 
rest of the particles are then assumed to form a kind of 
heat bath, which remains in equilibrium characterized 
by a temperature T = 1/(3 (which will be one of our 
parameters). All averages ((•■■)) are therefore thermal 
averages of the type: 

((...)) =Tr[(...)e-f" r \/Tr(e-f" T ). (4) 

However, we shall also allow for some external (to the 
above heat bath) influence acting on our system. There- 
fore we shall represent the operator b(k,t) as consisting 
of a part corresponding to the action of the heat bath, 
a(k,t), and also of a part describing action of these ex- 
ternal factors, R(k,t): 

b(k,t) = a(k,t)+R(k,t). (5) 



The time evolution of such a system is then assumed to 
be given by a Langevin equation [IQJ 

id t b(k,t) = F(k,t) - A(k,t) + P (6) 

(and a similar conjugate equation for b + {k,t)). These 
equations are supposed to model all aspects of the 
hadronization process (notice their similarity to equa- 
tions describing motion of Brown particle in some exter- 
nal field [Hj). The meaning of different terms appearing 
here is following: 

(i) The combination F(k,t) — A(k,t) represents the so 
called Langevin force and is therefore responsible for the 
internal dynamics of hadronization in the following man- 
ner: 

• A is related to stochastic dissipative forces and is 
given bv [Tnllflj 

/ + 0O 
drK(k,t-T)b{k,T), (7) 
-oo 

with the operator K(k,t) being a random evolu- 
tion field operator describing the random noise and 
satisfying the usual correlation-fluctuation relation 
for the Gaussian noise: 

(K + (k,t)K(k',t)) = 2^Wp~n8{k-k') (8) 

(k and p are parameters describing effect caused 
by this noise on the particle evolution in thermal 
environment jlflj). 

• The operator F(k, t) describes the influence of heat 
bath, 

F(k,t)= f +Q0 ^(k^c(k^e-^. (9) 

J-oo 27r 

(ii) Our heat bath is represented by an ensemble of 
damped oscillators, each described by operator c(fc M ) 
such that [c(fc M ), c + (fc^)] = <5 4 (fc M — k'), and character- 
ized by some function ^{k^), which is subjected to is a 
kind of normalization involving also dissipative forces 
represented by Fourier transformed operator K{k^), 
namely: 



r +co dcu 




Loo 27T 


K(k„)-LJ 



(iii) Finally, the constant term P (representing external 
source term in Langevin equation) denotes the possible 
influence of some external force (assumed here to be con- 
stant in time). This force would result, for example, in 
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a strong ordering of phases, leading therefore to the co- 
herence effect in the sense discussed in j2]. 

Out of many details (for which we refer to (§]) what 
is important in our case is the fact that the 2-particle 
correlation function for like-charge particles, as defined 
in (QJ, is given in such form = (u = k°,kj)): 



c 3 (Q) = OT- 

where 

/(fcp,^) = <fe + (fc M )6+(fcp6(fc M )6(fcp), 
/>„) = 



(11) 



(12) 

are the corresponding thermal statistical averages (in 
which temperature T enters as a parameter) with 
b(k tI ) — a(kfj,) + R(kfj,) being the corresponding Fourier 
transformed stationary solution of eq. 10 . As shown in 
P] (notice that operators i?(fc M ) by definition commute 
with themselves and with any other operator considered 
here) : 



+ (a+(k fl )a(k' li )}{a+(k^a(k fi )} + 

+ (a+ik^aik'^R+ik'JR^) + 

+ (h + (k'^(k^)R + (k^R(k^, (13) 

- (a+ikjai^)) + \R(k^)\ 2 . (14) 

This defines D(k fl , fc' ) in 11 111 in terms of the operators 
a(fc M ) and R(k^), which in our case are equal to: 



a(fc p 



and R(k^) 



(15) 



K(kJ-w K{k„)-u 

The multiplicity N depending factor £ is in our case 
equal to {(N) = (N) 2 /(N{N - 1)). This means there- 
fore that the correlation function Ci(Q), as defined by 
eq. ifTDl . is essentially given in terms of P and the two 
following thermal averages for the F(k,t) operators: 

{F+(k,t)F(k',t')) = 5 3 {k-k')- 

g |V| 2 n( W )e+-(*- t '), (16) 
{F{k,t)F + {k\t')) = S 3 (k-k')- 



(In; 
2^ 



H 2 [1 + n(w)]e- iw ( t -*') 



where = {exp [(u — [i)j3] — 1} 1 is the number of 

(by assumption - only bosonic in our case) damped os- 
cillators of energy ui in our reservoir characterized by 



parameters fi (chemical potential) and inverse tempera- 
ture (3 — l/T (both being free parameters). The origin of 
these parameters, the temperature = l/T and chemi- 
cal potential fi, is the Kubo-Martin-Schwinger condition 
that 

(a(k',t')a + {k,t)) = {a + {k,t)a{k' ,t-i/3))-exp(-P(x), (17) 

(see P1E3). This form of averages presented in llKit re- 
flects the corresponding averages for the c(fc^) operators, 
namely that 



(c(k,)c+(k' u )) 



S 4 (k tl -k')- n(ui) 



5\k 



*').[l + n(«)]. 



(18) 



Notice that with only delta functions present in 111 (it 
one would have a situation in which our hadronizing 
system would be described by some kind of white noise 
only. The integrals multiplying these delta functions and 
depending on (a) momentum characteristic of our heat 
bath ipikfj) (representing in our case, by definition, the 
hadronizing system) and (b) assumed bosonic statistics 
of produced secondaries resulting in factors n(u>) and 
1 + n(uS), respectively, bring the description of our sys- 
tem closer to reality. 

3. Results 

It is straightforward to realize that the existence of BEC, 
i.e., the fact that GziQ) > 1, is strictly connected with 
nonzero values of the thermal averages llllit . However, 
in the form presented there, they differ from zero only at 
one point, namely for Q = (i.e., for k^ = k'^). Actually, 
this is the price one pays for the QFT assumptions tac- 
itly made here, namely for the infinite spatial extension 
and for the uniformity of our reservoir. But we know 
from the experiment that C2{Q) reaches its maxi- 
mum at Q — and falls down towards its asymptotic 
value of C*2 = 1 at large of Q (actually already at Q ~ 1 
GeV/c). To reproduce the same behaviour by means of 
our approach here, one has to replace delta functions in 
eq. 11 fill by functions with supports larger than limited to 
a one point only. This means that such functions should 
not be infinite at — fc M — k'^ = but remain more 
or less sharply peaked at this point, otherwise remaining 
finite and falling to zero at small, but finite, values of 
\Qn\ (actually the same as those at which CziQ) reaches 
unity): 



Hk„ - =► ■ y/n(q = Q-r). 



(19) 



ISSN 0503-1265. Yk P . $ia. atcypn. 2003. T. 48, N 7 



3 



G.A.KOZLOV, O.V.UTYUZH AND G.WILK 



Here flo has the same dimension as the S function (actu- 
ally, it is nothing else but 4-dimensional volume restrict- 
ing the space-time region of particle production) and 
fl(q) is a dimensionless smearing function which con- 
tains the ^-dependence we shall be interested in here. 
In this way we are tacitly introducing a new parameter 
(mentioned already in the Introduction), r M , a 4- vector 
such that y/ {r^) 2 has dimension of length and which 
makes the product Q • r = Q^ry, = q dimensionless. This 
defines the region of nonvanishing density of oscillators 
c, which we shall identify with the space-time extensions 
of the hadronizing source. The expression l|19JI has to be 
understood in a symbolic sense, i.e., that fl(Q ■ r) is a 
function which in the limit of r — > oo becomes strict- 
ly a, 5 function. Making such replacement in eq. iflfil) 
one must also decide how to accordingly adjust n(uj) 
occurring there because now, in general, oj ^ w' . In 
what follows we shall simply replace n(u>) — > n{uj) with 
Oj = (lo+uj')/2 (which, for classical particles would mean 
that n(oj) — > y n(oj)n(cu')). 

In such way r becomes new (and from the phe- 
nomenological point of view also the most important) 
parameter entering here together with the whole func- 
tion fl(Q ■ r), to be deduced from comparison with ex- 
perimental data (one should notice at this point that the 
opposite line of reasoning has been used in ^2] where at 
first a kind of our Q(q) function was constructed for a 
finite source function and it was then demonstrated that 
in the limit of infinite, homogenous source one ends with 
a delta function) . With such a replacement one now has 



where 



fi(g) 



(1 + 00(1+0/) 



Cl(q) = r Sl(q), 7 = 



n 2 (u>) 
n(uj)n(uj') 



(20) 



,(21) 



with n{ui) the same as defined above. 

Another very important parameter entering 112011 is 
a. which first of all reflects action of external force P 
present in the evolution equation . This action is com- 
bined here (in a multiplicative way) with information 
on both the the momentum dependence of the reservoir 
(via \ip(k^)\ 2 ) and on the single particle distributions of 
the produced particles (via n(uj — /ixcoshy) where \it 
and y are, respectively, the transverse mass and rapid- 
ity). Parameter a summarizes therefore our knowledge 
of other than space-time characteristics of the hadroniz- 
ing source (given by Ct(q) introduced above). Notice that 



a > only when P^O. Actually, for a = one has 

K C 2 (Q) < l + 7fi(Q-r), (22) 

i.e., it is contained between limits corresponding to very 
large (lower limit) and very small (upper limit) values of 
P. Because of this a plays the role of the coherence pa- 
rameter 12] . For 7 ~ 1 , neglecting the possible energy- 
momentum dependence of a and assuming that a' = a 
one gets the expression 

c 2 {Q) = i + —^-^m + -^- I -m, (23) 

(l + o) z (l + o) z 

which is formally identical with eq. obtained in [2] by 
means of QS approach. It has precisely the same form, 
consisting two Q— dependent terms containing the infor- 
mation on the shape of the source, one being the square 
of the other, each multiplied by some combination of the 
chaoticity parameter p = 1/(1 + a) (however, in [2] p is 
defined as the ratio of the mean multiplicity of parti- 
cles produced by the so called chaotic component of the 
source to the mean total multiplicity, p = (JV c / l )/(JV)). 
In fact, because in general a ^ a' (due to the fact that 
oj ^ oj' and therefore the number of states, identified 
here with the number of particles with given energy, 
n(tu), are also different) one should rather use the gen- 
eral form lllll for C2 with details given by l!2(H and (12111 
and with a depending on such characteristics of the pro- 
duction process as temperature T and chemical potential 
fi occurring in definition of n{u>). 

Notice that eq. l12.'it differs from the usual empiri- 
cal parameterization of C2(Q) OP as given by eq.® in 
which < A < 1 is a free parameter adjusting the ob- 
served value of C2(Q = 0), which is customary called 
"incoherence "and with fl(Q ■ r) represented usually as 
Gaussian. Recently eq. I® has found strong theoretical 
support expressed in great detail in (where A is given 
by the same ratio od multiplicities as p above) . The nat- 
ural question arises: which of the two formulas presented 
here is correct? The answer we propose is: both are right 
in their own way. This is because each of them is based 
on different ways of defining coherence of the source. In 

one uses the notion of coherently and chaotically pro- 
duced particles or, in other words, one divides hadroniz- 
ing source into coherent and chaotic subsources. In [2] 
one introduces instead the notion of partially coherent 
fields representing produced particles, i.e., one has on- 
ly one source, which produces partially coherent fields. 
Our approach is similar as we describe our particle by 
operator b(k,t), which consists of two parts, cf. eq. ©, 
one of which depends on the external static force P. The 
action of this force is to order phases of particles in our 
source (represented by the heat bath). The strength of 
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this ordering depends on the value of the external force 
P. In any case, for P 7^ 0, it demonstrates itself as a 
partial coherence. 

Some comments are necessary at this point. Notice 
that it is of the same type as that considered in |2]- When 
comparing with [3 J one should notice that although our 
operators 3(fc /1 ) and i?(fc M ) look similar to operators de- 
fined in eqs. (4) and (5) of jH] they differ in the follow- 
ing. Our R{kfj) describes essentially the action of con- 
stant force P and as such it commutes with all other 
operators (including themselves). So it only introduces 
a partial ordering of phases of particles decreasing the 
C2 correlation function, i.e., acting as a coherent compo- 
nent, albeit we do not have coherent particles as such. It 
is also seen when realizing that in eq. l(T4*jl the two last 
terms contain only one pair of operators a. This in the 
language of [3 j translates to only one Wigner function, 
fch, to be present here. The operators R cannot form the 
second Wigner function (f co h in 0). This is the techni- 
cal origin of the three terms present in 1J2J (and in 
in comparison to two terms in and obtained in 0. 

Let us return to the problem of Q-dependence of 
BEC. One more remark is in order here. The problem 
with the S(k fl — k') function encountered in two particle 
distributions does not exist in the single particle distribu- 
tions, which are in our case given by eq. ([HI) and which 
can be written as /(fc M ) oc (a + (fc M )a(fc M )) + |i?(fc M )| 2 ~ 
(1 + a)(a + (/c M )a(fc /i )) (it is normalized to the mean mul- 
tiplicity: J d 4 kf(k li ) — (N)). To be more precise 



(24) 



where S(fc M , fc M ) is one-particle distribution function for 
the "free"(undistorted) operator a(fc M ) equal to 



n 



■0(fcp 



n(ui). 



(25) 



Notice that the actual shape of /(fc M ) is dictated both 
by n(u>) = n(cj;T, /x) (calculated for fixed temperature 
T and chemical potential fi at energy u as given by the 
Fourier transform of random field operator K and by 
shape of the reservoir in the momentum space provid- 
ed by ip(kfj,)) and by external force P in parameter a. 
They are both unknown, but because these details do 
not enter the BEC function C2(Q), we shall not pursue 
this problem further. What is important for us at the 
moment is that both the coherent and the incoherent 
part of the source have the same energy-momentum de- 
pendence (whereas in other approaches mentioned here 
they were usually assumed to be different). On the oth- 
er hand it is clear from that (N) = (N ch ) + (N coh ) 



(where (N c h) and (N co h) denote multiplicities of parti- 
cles produced chaotically and coherently, respectively) 
therefore justifying definition of chaoticity p mentioned 
above. 




a = 0.00 

a = 0.25 

a = 1.00 

a = 4.00 



0.0 0.2 0.4 0.6 0.8 1.0 



Q [GeV] 

1. Shapes of C2(Q) as given by eq. J2J - upper panel and for the 
truncated version of J2J (without the middle term) corresponding 
to eq. J3J - lower panel. Gaussian shape of £l(q) was used in both 



Fig. 1 shows in detail (using Gaussian shape of 0(g) 
function) the dependence of C^{Q) on different values of 
a = 0, 0.25, 1, 4 (again, used in the same approximate 
way as before and corresponding to p — 1., 0.8, 0.5, 0.2) 
and compare it to the case when the second term in eq. 
is neglected, as is the case in majority of phenomeno- 
logical fits to data. 



4. Summary and conclusions 

To summarize: using a specific version of QFT supple- 
mented by Langevin evolution equation © to describe 
hadronization process [HI E3 we have derived the usual 
BEC correlation function in the form explicitly showing: 



Figure 
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(i) - the origin of the so called coherence and its influ- 
ences on the structure of C2 (Q) ; (ii) - the origin of the (in- 
dependence of BEC represented by correlation function 
C2(Q). In our case the dynamical source of coherence is 
identified with the existence of a constant external term 
P in the Langevin equation. Its influence turns out to 
be identical with the one obtained before in the QS ap- 
proach 2\ and is described by eq. l|20)l. Its action is to 
order phases of the produced secondaries: for P — > 00 all 
phases are aligned in the same way and C2{Q) = 1- The 
coherence in j2] is thus property of fields and in our case 
property of operators describing produced particles. In 
both cases it occurs already on the level of a hadroniz- 
ing source. Dividing instead the hadronizing source itself 
into coherent and chaotic subsources results in eq. 
obtained in p|. The controversy between results given 
by j2] and pjj is therefore resolved: both approaches are 
right, one should only remember that they use different 
descriptions of the notion of coherence. It is therefore 
up to the experiment to decide which proposition is fol- 
lowed by nature: the simpler formula © or rather the 
more involved 11H together with 12011 . From Fig. 1 one 
can see that differences between both forms are clearly 
visible, especially for larger values of coherence a, i.e., 
for lower chaoticity parameter p. 

From our presentation it is also clear that the form 
of C2 reflects distributions of the space-time separation 
between the two observed particles rather than the dis- 
tribution of their separate production points (cf., for ex- 
ample, where it is advocated that it is in fact a 
Fourier transform of two-particle density profile of the 
hadronizing source, p(r\,r2) = p{r\ — r2), without ap- 
proximating it by the product of single-particle densities, 
as in 0). 

Finally, we would like to stress that our discussion 
is so far limited to only a single type of secondaries 
being produced. It is also aimed at a description of 
hadronization understood as kinetic freeze-out in some 
more detailed approaches. So far we were not interested 
in the other (highly model dependent) details of the 
particle production process. This is enough to obtain 
our general goals, i.e., to explain the possible dynamical 
origin of coherence in BEC and the origin of the specific 
shape of the correlation C'2(Q) functions as seen from the 
QFT perspective. We close with suggestion that both 
sources of coherence, that presented here and in [2] and 
that investigated in pj], should be considered together. 
The most general situation would be then hadronizing 
source composed with a number of subsources, each with 
different internal (discussed here) degree of coherence. 
For only one subsource present we would have situation 



described here whereas for a number of subsources, each 
being either totally coherent or totally chaotic, descrip- 
tion offered by [1,3, would automatically emerge. It is up 
to experimental data to decide and Fig. 2 tells us that 
it is not impossible task (at least in principle). 
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